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Abstract 

The Maxwell Garnett and Bruggeman formalisms were applied to estimate the effective permittivity 
dyadic of active dielectric composite materials. The active nature of the homogenized composite materials 
(HCMs) arises from one of the component materials which takes the form of InAs/GaAs quantum dots. 
Provided that the real parts of the permittivities of the component materials have the same sign, the 
Maxwell Garnett and Bruggeman formalisms give physically plausible estimates of the HCM permittivity 
dyadic that are in close agreement. However, if the real parts of the permittivities of the component 
materials have different signs then there are substantial differences between the Bruggeman and Maxwell 
Garnett estimates. Furthermore, these differences becomes enormous — with the Bruggeman estimate 
being physically implausible — as the imaginary parts of the permittivities of the component materials 
tend to zero. 

Keywords: Maxwell Garnett homogenization formalism; Bruggeman homogenization formalism; quantum 
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1 Introduction 

Suppose that we have a mixture of two isotropic dielectric materials, specified by the relative permittivity 
scalars t\ and 62. Provided that the length scale of nonhomogcncitics in the mixture is small compared with 
electromagnetic wavelengths in both component materials [1] , the mixture can be treated as a homogenized 
composite material (HCM). The familiar formalisms named in honour of Maxwell Garnett [2, 3, 4] and 
Bruggeman [5, 6] provide a means of estimating the relative permittivity of the HCM. 3 These formalisms are 
well-established for dissipative HCMs arising from component materials for which the real parts of ei and 
62 have the same sign [8]. However, we recently demonstrated that these formalisms give rise to estimates 
of the HCM's relative permittivity which are not physically plausible if the component materials are weakly 
dissipative and the real parts of ei and £2 have different signs [9, 10]. Furthermore, this limitation also 
extends to the Bergman-Milton bounds on the HCM's relative permittivity [11]. 

In this short communication, we consider whether these limitations on conventional homogenization 
formalisms also extend to active HCMs. This topic is particularly timely, given the current huge interest 
in metamaterials — many of which can be viewed as HCMs [12]. While certain metamaterials can be 
engineered to exhibit negative refraction, the issue of dissipation presents a major barrier to their practical 
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implementation. In order to overcome inherent dissipation, negatively refracting metamaterials containing 
active components have been proposed recently [13, 14]. 

In the following, double underlining signifies a 3x3 dyadic, with and 7 denoting the null and identity 
dyadics respectively. The unit vector aligned with the Cartesian z axis is written as z. The real and imaginary 
parts of complex quantities are identified by the prefixes 'Re' and 'Im', respectively; and i = >/— T. 



2 Investigation 

2.1 Component materials 

We base our investigation on the metal-wire/quantum-dot composite material proposed by Bratkovsky et 
al. [14] . That is, we consider a three-component composite material comprising silver nanorods inserted into 
an array of InAs / GaAs quantum dots, with the assembly being immersed in a host dielectric material. The 
silver nanorods are aligned with the Cartesian z axis, while the quantum dots are presumed to be spherical 
particles. The particles which comprise the host dielectric material are assigned a spherical shape in the 
implementation of the Bruggeman homogenization formalism. In contrast, no microstructure is assigned to 
the host dielectric medium in the implementation of the Maxwell Garnett homogenization formalism. The 
volume fractions of the silver nanorods, quantum dots and host dielectric material arc written as f^ g , /qd 
and fh, respectively, and we have 

fAg + fQD + fh = 1. (1) 

The relative permittivity of the silver nanorods is expressed as a function of angular frequency u) by the 
Drude formula [15] 

eA fl H = l--2-^ , (2) 

where the plasma frequency u p = 1.38 x 10 16 rad s _1 and the relaxation rate jA g = 10 14 s _1 . The relative 
permittivity of the quantum dots is provided as [14] 

£QDM = efcH • (3) 

uj - Ll>„ + l "f QD 

Herein, hw = 0.8 eV, ha = 22.9 meV, where H is the reduced Planck constant, the background relative 
permittivity ej, = 11.8, and the broadening parameter jqd = a/15.3 [16]. The host dielectric material is 
taken to have the frequency-independent relative permittivity e/j = 3. The real and imaginary parts of 6A g 
and €qd are plotted against free-space wavelength A Q e [1, 2] /im in Fig. 1. We see that lu\eq D < which 
indicates the quantum dots are active entities, particularly so at the resonance centred on A Q = 1.55 /im. 
From the point of view of the applicability of homogenization formalisms, there are two key points to notice: 
(i) the quantity 

S= ^L (4) 
Rce QD 

is negative over the entire wavelength range considered; and (ii) we have | Ree^g | S> | ImtA 9 | for A D e [1, 2] 
/im, and | ReeQD | 3> | Imegc | for A Q G [1, 1.55 — rj] U [1.55 + 77, 2] /im where 77 specifies a sharp resonance 
region centred on A Q = 1.55 /im. By analogy with dissipative HCMs [9, 10, 11], points (i) and (ii) suggest 
that the application of homogenization formalisms may be problematic. 

2.2 Homogenization formalisms 

The HCM is characterized by the relative permittivity dyadic 

g HCM = e HCM ( £ _ £ z) + £ H CM ~ ~_ (5) 

We write 'MC or 'Br' in lieu of 'HCM' according to whether the Maxwell Garnett or the Bruggeman estimate 
of e HCM is being considered. The anisotropy of the HCM arises due to the orientation of silver nanorods. 
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The Maxwell Garnett estimate of the relative permittivity dyadic is given explicitly by [17] 



| MG = e h L + f Ag g A9/h ' ( L ~ ^ * Aalh 



wherein the polarizability density dyadics 



f QD g QD/h • [L-^-g QD/h 



3e/» - 



(£ = Ag,QD), 



> . 



i jh = (ee - e h ) [L+ (et - e h ) R e/h 
and the depolarization dyadics [18] 

D Ag/h = J_ tj_ ii \ 

= 2e h V = ; 

D QD/h = 1 j 

= 3e/j = 

The Bruggeman estimate of the relative permittivity dyadic is given implicitly by [17] 

f Ag g Ag/Br + fan g QD/Br + h g h/Br = £ , 
wherein the polarizability density dyadics 

= ( e ,| _ ^r) . ^ + g/Br . ^ _ 1 {£ = ^ QD ^ 

and the depolarization dyadics [18, 19] 

D A ^ Br = -L- (1-zz) 

E t/Br = -BT [aL(a)(L-zz)+L z (a)zz] , (l = QD,h) 
Herein, the scalar functions 

L{a) = ^F{a)-^— 1 [l-F{a)\ 

L z (a) = -^—[l-F(a)] 

a — 1 



with 



and the scalar a = ef r /ef r . 



F(a) = . tan" 1 Wa^T) 
va — 1 



(6) 
(7) 

(8) 

(9) 
(10) 

(11) 



(12) 



(13) 



The relative permittivity dyadic e Br may be extracted from the nonlinear equation (9) by numerically. 
For example, the iterative scheme [7, 20] 



e Br =T (Br 1 

= 0) l=0-i)J 
wherein the operator T is defined as 



U = 1.2,3,. 



(14) 



T{g Br } - {fA g e Ag [L + R A3/Br ' {eA g L-g Br )Y 1 +fQDe QD [L + R QD/Br ' (e Q DL-g 

-l 



+f h e h [l = + a h ' Br • (e h L-g Br )Y } ' {/a 9 [L + R A9/Br ' (e Ag L-f 
fan [l + R QD/Br ' (cQnl-g^)]' 1 +fh[L + R h/Br ' (^I-i^rT 



and the initial value 



= B (Q) = ( -^f 6 ^ + fao^QD + fh£h ) I , 



(15) 



(16) 



works sufficiently well for our purposes here. 
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2.3 Numerical estimates of the HCM permittivity 



Let us begin by considering the Maxwell Garnett and Bruggeman estimates for an active HCM when the 
parameter 5 is positive. This is achieved by setting f Ag = 0. Thus, the composite material comprises only 
quantum dots and the host dielectric material, and the resulting HCM is an isotropic dielectric material 
with relative permittivity dyadic e HCM L Estimates of the HCM relative permittivity computed using the 
Maxwell Garnett formalism and the Bruggeman formalism, namely e MG and e Br respectively, are plotted 
against A Q in Fig. 2, when Jq D = 0.3. The Maxwell Garnett and Bruggeman estimates are similar, both 
qualitatively and quantitatively. For both formalisms, the HCM is active across the entire wavelength range 
considered, and its permittivity scalar exhibits a sharp resonance at A Q = 1.55 yum. Therefore, we infer that 
both the Maxwell Garnett and the Bruggeman formalisms are suitable for active HCMs with 5 > 0. 

Next we turn to the 5 < regime which is known to be problematic for dissipative HCMs [9, 10, 11]. In 
keeping with Bratkovsky et al. [14], we fix fAg = 0.063 and /qd = 0.3; thus, fh is slightly reduced from 
its value for Fig. 2. In Fig. 3, the two entries in of the HCM's relative permittivity dyadic, as computed 
using the Maxwell Garnett and Bruggeman formalisms, are plotted against A . While the Maxwell Garnett 
and Bruggeman graphs in Fig. 3 are qualitatively similar, there are substantial quantitative differences. 
In particular, the imaginary part of ef r is more than 20 times larger than e^ G , apart from at the small 
resonance region centred on A Q = 1.55 fxm. There are also substantial differences between Re ef r and Re 
ef G - By comparison, the differences between ef r and e^ G are relatively small. 

For dissipative HCMs, the problematic nature of the 6 < regime is exacerbated as the imaginary parts 
of the constitutive parameters tend towards zero [9, 10, 11]. In order to investigate this issue for active 
HCMs, we repeated the computations of Fig. 3 with jAg and jqd multiplied by 10 -6 . The corresponding 
plots of the real and imaginary parts of cfl' M are presented in Fig. 4. The differences between the Maxwell 
Garnett estimates and the Bruggeman estimates are now enormous, especially between the imaginary parts 
of ef- andef G . 



3 Concluding remarks 

Our numerical studies reveal that the well-known homogenization formalisms named after Maxwell Garnett 
and Bruggeman appear to be fairly consistent, and provide physically plausible estimates, when applied to 
active isotropic dielectric HCMs, provided that the real parts of the relative permittivities of the component 
materials have the same sign. In contrast, if the real parts of the permittivities of the component mate- 
rials have different signs then there are substantial differences between the estimates yielded by the two 
formalisms. These differences become enormous in the limiting case in which the imaginary parts of the 
relative permittivities of the component materials become vanishingly small. 

For the particular homogenization scenario considered here, in the 5 < with jAg,QD — * regime, 
the HCM permittivity parameters estimated using the Bruggeman formalism have relatively large positive- 
valued imaginary parts away from the resonance region centred on A Q = 1.55 /xm, and relatively large 
negative-valued imaginary parts at the resonance region centred on A Q = 1.55 /xm. This is not physically 
plausible since it implies that the HCM is a strongly active material at the resonance region, and a strongly 
dissipative material away from the resonance region, in the limit wherein the component materials become 
inactive and nondissipativc. Let us note that Bruggeman formalism arises as the lowest-order formulation 
of the strong-pcrmittivity-fluctuation theory (SPFT) 4 [22]. Accordingly, higher-order implementations of 
the SPFT are also subject to the limitations of the Bruggeman formalism. Furthermore, since the Maxwell 
Garnett and Bruggeman formalisms share a common provenance [23] , doubt is also cast over the applicability 
of the Maxwell Garnett formalism for active HCMs with d < where the imaginary parts of the relative 
permittivities of the component materials are relatively small. 

The present study extends and reinforces our previous studies which highlighted the limitations of con- 
ventional homogenization formalisms when applied to scenarios wherein the real parts of the constitutive 

4 Otherwise known as the strong-property-fluctuation theory for more general HCMs [21]. 
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parameters characterizing the component materials have different signs [9, 10, 11]. We see here that caution 
is needed for active HCMs as well as dissipative ones. 



Acknowledgement: TGM thanks Dr Petter Holmstrom (KTH Royal Institute of Technology, Sweden) for 
his helpful comments. 
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Figure 1: Real and imaginary parts of the relative permittivities tAg and £qd plotted against free-space 
wavelength A c (pm). 



(i 



J 


Ag = o 


\ 


/qd = 0.3 


1.0 1.2 1.4 1.6 1.8 2.0 
l (/mi) 




) - 
) - 


/qd = 0.3 


1.0 1.2 1.4 


1.6 1.8 2.0 



a 4.5 



J 


/Ag = 




/qd = 0.3 



1.0 1.2 1.4 1.6 1.8 2.0 

A„ (/jm) 




Figure 2: Bruggeman (left) and Maxwell Garnett (right) estimates of the HCM's relative permittivity plotted 
against A (/im), with fA g — and Jqd = 0.3. 
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Figure 3: Bruggeman (left) and Maxwell Garnett (right) estimates of the two entries in the HCM's relative 
permittivity dyadic plotted against A Q (/im), with f^ g — 0.063 and Jqd = 0.3. The solid curves represent 



r, whereas the dashed curves represent ef r,MG . 
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Figure 4: As Fig. 3 but with both ^Ag and ^qd multiplied by a factor of 10 
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